Quantum damping of Fermi-Pasta-Ulam revivals in ultracold Bose gases 
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We propose an experimental scheme for studying the Fermi-Pasta-Ulam (FPU) phenomenon in a 
quantum mechanical regime using ultracold atoms. Specifically, we suggest and analyze a setup of 
one-dimensional Bose gases confined into an optical lattice. The strength of quantum fluctuations 
is controlled by tuning the number of atoms per lattice sites (filling factor). By simulating the 
real-time dynamics of the Bose-Hubbard model by means of the exact numerical method of time- 
evolving block decimation, we investigate the effects of quantum fluctuations on the FPU recurrence 
and show that strong quantum fluctuations cause significant damping of the FPU oscillation. 

PACS numbers: 03.75.Kk, 03.75.Lm 
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Systems of ultracold atomic gases in optical lattices 
have offered unique possibilities to study nonequilibrium 
dynamics of quantum many-body systems [if Q . Ex- 
treme cleanness and exquisite controllability of cold atom 
systems have led to experimental realization of various in- 
teresting nonequilibrium phenomena, such as the dynam- 
ics following the quench across the supcrfluid (SF)-Mott 
insulator (MI) transition Q, the significant suppression 
of transport of onc-dimensional (ID) lattice bosons [3, HI, 
and the nonergodic dynamics of (nearly-)integrable ID 
Bose gases Q. In particular, the experiment of Kinoshita 
et al. [6( has shed light on the fundamental problem of the 
statistical mechanics, which is the relation between inte- 
grability and ergodicity, and stimulated renewed theoret- 
ical interest in this topic (see Ref. Q and Refs. therein). 

It is widely known through the celebrated numerical 
experiment by Fermi, Pasta, and Ulam (FPU) that non- 
integrable systems can also exhibit nonergodic dynamics 
under a certain condition 0. In order to address the 
essential question whether the nonlinearity leads to er- 
godicity, FPU studied dynamics of ID chains of classical 
nonlinear oscillators where initially only the first normal 
mode is excited. Contrary to their expectations, FPU 
found quasi-periodic recurrences of the initially excited 
mode, which are now called FPU recurrences. This dis- 
covery is marked as one of the most important milestones 
in the nonlinear physics because it largely stimulated the 
whole development of the chaos theory Q. There have 
been multiple attempts to experimentally demonstrate 
FPU recurrences in several systems, including electrical 
networks Q , double-plasma devices [l(| , deep water [Tljj , 



magnetic films [12|, [13(, and optical fibers [14|. How- 
ever, such experiments usually suffer from the difficulty 
in making a physical system isolated from sources of en- 
ergy dissipation that preclude the demonstration of com- 
plete recurrence 0-12|. 

Theoretically, there does not appear to be a clear rea- 
son why FPU phenomenon, i.e. revivals of non-integrable 



many-body dynamics, should not be a general feature 
in low dimensional non-equilibrium physics, not particu- 
larly confined to the model studied by FPU. Indeed, in 
this Letter, we propose an experimental scheme for ob- 
serving the FPU recurrences with ultracold Bose gases 
in optical lattices and present the numerical analysis of 
the corresponding dynamics. Ultracold atomic systems 
are especially suited for thorough studies of the FPU dy- 
namics because of their unprecedented controllability and 
isolation. The central idea to achieve this goal is two- 
fold. First, when the filling factor is sufficiently large, the 
Bose-Hubbard (BH) model describing ultracold bosons 
in an optical lattice exhibits features of the FPU model, 
such as recurrences and energy localization in q-space. 
Secondly, in this regime, one can generate the initial 
state with the first normal mode excited prominently by 
means of phase imprinting. We show that dynamics of 
f D Bose gases in optical lattices exhibit the FPU recur- 
rences in the classical (Gross-Pitaevskii) limit. Moreover, 
such setup has another advantage allowing one to control 
the strength of quantum fluctuations through varying the 
filling factor, the interatomic interaction, and the lattice 
depth. To our knowledge this is the first experimentally 
feasible proposal to study the FPU phenomenon in the 
quantum regime. We also analyze dynamics in the sug- 
gested setup usin g t he time-evolving block decimation 
(TEBD) method We find that the quantum fluc- 

tuations cause damping of the FPU recurrences. Such 
damping can be reduced by increasing the filling factor 
such that in the classical limit we find nearly complete 
revivals characteristic for the original FPU problem. 

We consider a system of N bosons in a ID optical 
lattice. Assuming that the lattice is sufficiently deep, the 
system can be described by the ID BH Hamiltonian [lij : 
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where the field operator St (Sj) creates (annihilates) a 
boson on the j-th site, J is the hopping energy, and 
U is the onsite interaction. To make closer connection 
with the FPU's original model, we assume the absence 
of a harmonic trap that is used in standard experiments. 
This assumption can be valid when one uses a box-shaped 
trap 17 1 or a ring optical lattice [HI . We henceforth fo- 
cus on the case of a ring geometry for theoretical simplic- 
ity taking the periodic boundary conditions in Eq. (TTJ. 

The strength of quantum fluctuations is characterized 
by the effective Planck's constant h e = y/U/(yJ), where 
v = N/L is the filling factor. Physically h c controls prox- 
imity of the system to the Mott insulator regime with the 
transition occurring at h c ~ 1.5 for integer lar ge v . h e 
can be experimentally varied in a broad range 1|1, 20j [. 
We note that the nonlinearity in the system is controlled 
by another parameter k = Uv/J. When this parame- 
ter is large, the number fluctuations in the system are 
small and the Hamiltonian (TTJ can be mapped into the 
quantum rotor model (l9| . 

Let us first consider the classical, Gross-Pitaevskii, 
limit (h e — > 0) and identify parameter regions where FPU 
recurrences occur. In this limit, one can safely neglect 
quantum fluctuations and replace the field operator bj (t) 
with a classical field j (t) . Carrying out the replacement 
in Eq. (TTJ), we obtain the classical Hamiltonian: 
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FIG. 1: (color online) (a) The energy of all the normal modes 
at t = 0, where A = 0.4. (b)-(d)Time evolution of the energy 
of the first four modes A = 0.05 (b), 0.4 (c), and 2.0 (d). We 
set L — 32 and n = 20. 



The dynamics of the system can be described by the dis- 
crete nonlinear Schrodinger equation (DNLSE), or equiv- 
alently discrete Gross-Pitaevskii equation [2lj |: 



ih — - = ~~z — = -J(Vj- 
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It is clear that in this classical limit both static and dy- 
namic properties of the DNLSE (after trivial time rescal- 
ing by J) are determined by the parameter k = Uv/J 
introduced earlier. In the rotor regime k S> 1, DNLSE 
can be mapped into 



dt 2 
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present in Eq. (j4j which is equivalent to an additional 
small cubic nonlinearity in the FPU model. Indeed, we 
work in the regime where cubic nonlinearity is not neg- 
ligible. We stress that this is perfectly fine for studying 
FPU phenomena since the only necessary condition is 
for k to be sufficiently large, and we do not need to be 
strictly in the rotor limit. 

In order to induce FPU recurrence, we consider the 
following initial condition, 



4>j(t = 0) = A sin (2-rrj/L) , nj(t = 0)=u, 



(5) 



which can be g enerated with the use of phase-imprinting 
techniques [22| . The amplitude A is related to the energy 



where ^ and n 3 - are the phase and the density of the con- E A = H(t = 0) - U gs injected by this phase imprint as 



e^i and <fj t i = 4>j — <j>{. No- 



densate defined by \&j = 
tice that the condition k 3> 1 is compatible with h c -C 1 
at large filling factors i/> 1. DNLSE in the rotor regime 
reduces to the quartic FPU model when the nonlinearity 
is small: \<pj — 4>j+i\ < 1, so that one can expand the 
sinusoidal terms as sin ip ~ ip— ip 3 /6. However, there are 
also some differences; for example, the DNLSE conserves 
the number of particles, ^ \^j\ 2 = N, while there is no 
such conservation in the latter models. In addition at 
large but finite n there is a small quadratic nonlinearity 



Ea — 2NJ(ttA/L) 2 , where H gs is the classical ground 
state energy. Fig. [TJa) shows the initial energy distribu- 
tion of the normal modes E q (t = 0) for the chosen initial 
condition <j5j) for L = 32 and K = 20. The normal mode 
energies are given by 



E q = 4isJs'm- { —q 
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In Eq. ©, 

of the phase 
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are the Fourier transforms 



and the density n, given by 
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FIG. 2: (color online) Time evolution of the momentum occu- 
pation n q (t) for q — 1 (blue solid line). The parameters of the 
system are the same as in Fig. [He). The red dashed line rep- 
resents C x E q -i(t), where the constant C is chosen to agree 
with the envelope of the oscillation of n q= i(t). The graphs 
(a) and (b) show short and long time dynamics respectively. 



L- 1 / 2 J2j<f>je~ i2nqj/L ,n q = L-V' 2 Y,j n 3 eri2 ' K ' l5/L i and 
the mode number q is integer satisfying 1 < q < L — 1. 
In Fig. [TJa) , it is clear that only the normal modes with 
the lowest frequency (q = 1 and L — 1) are prominently 
excited in the initial state as in the original work of FPU. 
Therefore, the FPU recurrences are expected to occur in 
the dynamics subjected to the initial phase imprint. 

In Figs.QJb)-(d), we show time evolution of the energy 
of the first four modes for several values of the ampli- 
tude. These results are obtained by numerically solv- 
ing DNLSE ([3]). When the amplitude A is small as in 
Fig. |Tfb) , the energy of the initially excited mode is al- 
most constant because the coupling among the normal 
modes is too weak to induce the mode mixing needed for 
the FPU dynamics. As A increases, the energy spreads 
to other modes but still remains confined to a few of 
them. In this regime one observes repetitive recurrences 
of the first mode (Fig. Die)). Both long-time revivals 
and even longer time super-revivals, where nearly 100% 
of the energy returns to the first mode, characteristic for 
the FPU dynamics Q are clearly present. Finally as A 
increases even further and exceeds a certain threshold, 
the dynamics becomes chaotic rapidly reaching equipar- 
tition of the energy between all the modes (Fig. Hfd)). 
Such a threshold for ergodic behavior has been previously 
found in many other models including the quartic FPU 
model [23, 24 1, the nonlinear Schrodinger equation with 
open boundaries [25| , and DNLSE at smaller filling [2(| . 

The temporal behavior of the mode energies shows 
close analogy between BH and FPU models. However 
these mode energies are not easily observable in cold 
atom experiments. It is thus much more natural to ana- 
lyze time evolution of the quasi-momentum occupation: 



= L 



3,1 *j*' e " 



= l27T iU-i)/ L i which is a standard ob- 
servable [1[. l n fig-HI a ) we plot n q —i(t). Clearly the mo- 
mentum distribution oscillates at much larger frequency 



FIG. 3: (color online) Time evolution of the energy E q of the 
first four modes, where L — 8 and Uv/J — 20. (a) The clas- 
sical dynamics based on DNLSE. (b) The quantum dynamics 
calculated by applying TEBD to the Bose-Hubbard model. 
We set v = 10, which means that the effective Planck's con- 
stant h e = \/0T2 ~ 0.447. 



■if 



u w \/2vJlJ-n /(Lh) compared to frequency of the en- 
ergy oscillations. These oscillations are expected even 
in the absence of the coupling between normal modes. 
From fig. [2Kb) it is obvious that the envelope of these 
fast oscillations exhibits the FPU recurrences and thus 
the quasi-momentum distribution can serve as a good 
observable for detecting the FPU dynamics. 

Having confirmed that the FPU recurrences occur be- 
low the equilibration threshold in the classical limit, we 
now consider quantum dynamics of our system. For 
this purpose we apply the quasi-exact numerical method 
of TEBD [l5| with periodic boundary conditions [27J 
to the BH model. For efficient TEBD simulations for 
large filling factor [y 3> 1) we reduce the dimension of 
the local Hibert space by introducing a lower bound of 
the particle occupation number in addition to an up- 
per bound [2(J. An initial state for the quantum FPU 
dynamics is generated by preparing the system in the 
ground state of the Hamiltonian ([1]) via the imaginary 
time propagation and then applying the phase imprint 

operator, P = Y[j=i e~ l8jb i bj ', to the ground state at 
t = 0, where 0j = Asm{2irj / L), As in the classical 
case, we calculate the time evolution of the mode ener- 
gies E q (t) additionally subtracting the contribution from 
the depletion which is equal to the energy of the 
corresponding mode in the ground state. 

In Fig. [31 we show time evolution of the mode energies 
both for the classical (a) and the quantum (b) dynam- 
ics. In the quantum case we use v = 10 corresponding 
to h e ~ 0.447. In both cases n — 20 and L = 8. Note 
that L = 8 is a minimum possible size required to ob- 
serve the FPU dynamics. For smaller systems the fre- 
quency of the FPU oscillation is too close to that of the 
lowest-energy excitation. For even this relatively small 
system the Hilbcrt space is too large to allow for ex- 
act diagonalization. From Fig [3jL) it is clear that the 
FPU oscillations are significantly damped in the quan- 
tum case. To corroborate that this damping is due to 
quantum fluctuations, we vary the filling factor v while 



4 




10 15 
tJ/% 



FIG. 4: (color online) Time evolution of the energy E q -\ of 
the first mode for v = oo (DNLSE), 20, 14, 10, and 8 where 
L = 8 and Uu/ J — 20. The inset depicts the damping rate V 
of the FPU oscillation as a function of the effective Planck's 
constant h c = ^JU/JyJ). 



fixing k = 20, and analyze how the damping depends 
on the effective Planck's constant h c . In Fig. 2] we show 
the time evolution of the first mode energies for several 
values of v. Clearly the damping becomes stronger at 
smaller v, i.e. larger h c . We extract the damping rate Y 
by fitting E q= \{t) to a function f(t) = ae~ Tt cos(fei) + c, 
where a, b, c, and Y are free parameters. In the inset 
of Fig. |4l we plot the extracted damping rate versus h c . 
When the filling v becomes even smaller so that the sys- 
tem approaches the Mott insulator phase we find that 
the mode energies start exhibiting new quantum revivals 
which are not related to the FPU physics, but rather to 
the gradual emergence of the particle-hole excitations. 

It is worth noting that the damping of the FPU oscilla- 
tion does not result in equilibration within the timescale 
of our TEBD calculations. This implies that quantum 
fluctuations do not necessarily make the FPU dynamics 
crgodic. Instead, the system seems to approach a certain 
metastable state, and it is interesting to seek a statis- 
tical mechanics describing such a metastable state. We 
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will address this issue in detail elsewhere 

One possible qualitative explanation for the damping 
can be made from the perspective of the decay of mater- 
wave solitons due to quantum fluctuations. On the basis 
of the soliton concept, the classical FPU recurrences are 
interpreted as follows [29|]. The initially prepared wave 
with the sine form splits into a number of solitons with 
different velocities. Since each soliton is stable during 
time evolution, the solitons come back very close to the 
initial configuration after colliding many times, produc- 
ing the recurrences. In the quantum case, on the con- 
trary, it was argued that the solitons arc unstable 
and that the soliton lifetime decreases with strengthen- 
ing quantum fluctuations [loj]. This soliton decay may 
be responsible for the damping of the FPU oscillations. 

In conclusion, we analyzed the dynamics of one dimen- 
sional bosons in optical lattices and demonstrated that 
under initial phase imprint at high filling factors FPU re- 
currences occur in this system. We numerically studied 



time evolution of this system to reveal quantum effects 
on FPU recurrences. We showed that the damping of the 
FPU oscillations becomes more significant with increas- 
ing the strength of quantum fluctuations. We emphasize 
that the controllability of quantum fluctuations in ultra- 
cold atom systems can open up a new direction of the 
study of the FPU phenomenon. 
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